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Abstract. We present a theoretical description of a system of many 
spins strongly coupled to a bosonic chain. We rely on the use of a spin- 
wave theory describing the Gaussian fluctuations around the mean-field 
solution, and focus on spin-boson chains arising as a generalization of 
the Dicke Hamiltonian. Our model is motivated by experimental setups 
such as trapped ions, or atoms/qubits coupled to cavity arrays. This 
situation corresponds to the cooperative (Eg)/3) Jahn- Teller distortion 
studied in solid-state physics. However, the ability to tune the parame- 
ters of the model in quantum optical setups opens up a variety of novel 
intriguing situations. The main focus of this paper is to review the 
spin-wave theoretical description of this problem as well as to test the 
validity of mean-field theory. Our main result is that deviations from 
mean-field effects are determined by the interplay between magnetic 
order and mesoscopic cooperativity effects, being the latter strongly 
size-dependent . 



1 Introduction 

Experimental progress in quantum optical setups has opened up new research ar- 
eas where the controllability of those experimental systems meets the complexity of 
quantum many body physics. For example, trapped ions and ultracold atoms can 
now be controlled to the extent in which they emulate the physics of Condensed Mat- 
ter systems in a way that may lead to the implementation of quantum simulators 
of many-body models |l|2j. The latter are devices in which quantum states may be 
prepared and measured and interactions may be tuned, so as to mimic complex dy- 
namics with a practical and scientific interest in material science. The research field 
of Analogical Quantum Simulation has emerged to exploit this idea. 

Remarkably the physical elements of many of those systems, particularly trapped 
ions and atoms interacting with photons, can be understood by means of spin-boson 
models. Effective spins are implemented by two-level systems which can be either 
atomic transitions of trapped ions |3l4j or isolated energy levels in solid-state qubits 
[617] . The bosonic degree of freedom is provided by vibrations in ion Coulomb crystals 
[8] or photons in optical or microwave cavities. Bosonic modes may be collective 
excitations either because of the vibrational couplings in the case of ions, or because 
of the collective nature of photonic bands in cavities. Finally, a number of different 
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spin-boson couplings may be implemented in those systems, either by laser-induced 
forces in the case ions, or by the light-matter coupling between qubits and cavities. To 
focus our discussion, let us define crj — — |0)j(0|,with |0) 3 , 1 1)^- internal energy 

levels of an atom (qubit), and ctj, annihilation, creation operators for a bosonic 
mode at site j. A natural coupling in quantum optical systems is 

3 

If an additional process couples different bosonic modes, qubit-boson couplings pro- 
vide us with a mechanism to couple different qubits. This idea is the basis for quantum 
gate designs |9I7I10) . In the field of Analogical Quantum Simulation, the use of cou- 
plings like Q has been mainly an instrument to get effective models, like quantum 
magnetic Hamiltonians, in which tracing out bosons results in effective interactions 
of the form J2j k ^fe <7 l cr I [415111112113114) . Also, strongly coupled spin-boson chains 
have been considered in a few previous works, for example in the context of disor- 
dered systems |15j and quantum dissipation involving a single spin interacting with 
a bosonic bath [16 . 

Quite remarkably, however, interaction (fT]) poses an outstanding many-body model 
defined on a quantum lattice, with a quantum phase diagram which is intrinsically 
interesting. It belongs to a family of many-body models previously considered in 
the context of solid-state physics, in particular, in the description of the coopera- 
tive Jahn- Teller E<g>/3 distortion [17] . The latter is a many-body effect which appears 
in solids, where two-level systems correspond to orbital electronic levels localized in 
atoms, which interact with the lattice vibrations. This effect has received attention 
for its connections with high-Tc superconductivity and colossal magnetoresistance 
[181191201 . Very recently it was shown that a generalization of interaction (JTJ) can 
be implemented in a quite natural way in a linear ion crystal placed in a magnetic 
field gradient [21122) . The prospects for experimental investigation of the cooperative 
Jahn- Teller effect would be very interesting, since they would allow us to observe 
novel quantum structural phase transitions in clean controllable systems. Solid-state 
systems pose many limitations to observe quantum effects in cooperative Jahn- Teller 
models, since experiments are typically performed at high temperatures, and also 
because Jahn- Teller couplings are masked by magnetic or electronic dynamics. A va- 
riety of related models are currently under theoretical investigation by several groups 
including the single particle case [23], Jaynes-Cummings- Hubbard models [24125] . 
structural phase transitions in quantum potentials [23] and spin-Peierls transitions 

This paper is dedicated to the spin- wave theory of cooperative (E<g)/3) Jahn- Teller 
models arising from longitudinal a z couplings of the form (flj. A natural starting 
point is mean-field theory together with a spin-wave description of Gaussian fluctu- 
ations [25]. The main goal of this paper is to test the validity of mean- field theory 
by calculating the intensity of quantum fluctuations. The outline of the paper is the 
following. In section 2 we introduce the Hamiltonian that we are considering for co- 
operative Jahn- Teller systems, and present some qualitative considerations arising 
from the comparison with the Dicke model. In section 3 we review the mean-field 
in this model. In section 4 we present a calculation of the Gaussian fluctuations by 
means of the definition of Holstein-Primakoff bosons. In section 5 we study in detail 
the case of Periodic Boundary Conditions and check the validity of the mean-field 
approximation in a region that depends both on the magnetic order of the system, 
as well as on mesoscopic cooperative effects as a function of the number of particles. 
Finally in section 6 we present our conclusions. 
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2 Cooperative Jahn-Teller systems 

We consider a chain of N spins coupled to N bosonic modes. In addition, there are 
couplings between different bosonic modes as well as a magnetic field acting on the 
spins. This leads to a cooperative Jahn-Teller Hamiltonian of the form, 

Hjt = H s + Hb + H s b, 

3 3 3,1 3 

where aj, aj are annihilation/creation operators of bosons that are localized near the 
effective spin. The physical meaning of those bosonic modes can be either vibrations 
of a trapped ion, or photons confined in a cavity, fl is a transverse magnetic field, and 
ojj and tjj are boson local energies and hopping amplitudes; the latter determine the 
boson dispersion. There are several experimentally relevant cases. However, in this 
work we focus on nearest-neighbours boson tunneling, 

= -*(<W + S 3,i-i), * > °- ( 3 ) 

The latter describes coupled cavities, and to a fair approximation also the radial 
phonons of an ion Coulomb chain |29j . We diagonalize the bosonic Hamiltonian by 
transforming to a collective mode basis, a,j — ^2 n Mj_ n a n , where Mj <n is the collective 
mode amplitude at site j. The relation 

U n 5n,m = X! MtnitjJ + u j^j,l) M l,m (4) 
3,1 

determines the collective mode energies, Q n , such that = ^ n uina^cin. Finally, the 
spin-boson coupling in the collective modes basis reads 

H sh - g J2 °3 ( M ^a n + Ml n ai) . (5) 

The problem posed by £Tjt resembles the Ising quantum magnet in a transverse 
field. For large fl the ground state is a product state of the boson vacuum and spins 
pointing in x. In the case of large g magnetic order is in the ^-direction. To show this 
result, let us see how to get effective spin-spin Ising interactions |4I5] . Consider the 
canonical transformation U = e~ s with 

S = Y^g/u^o*} (M j>n a„ - M* n al) , (6) 
representing a spin-dependent boson displacement. We apply U to the fl = case, 

„2 

e- s (H h + H sh )e s = H h + H sh = Y, w n o£a„ - ^ ^ M tn^l M l,n, (7) 

n j,l, n U ' n 

which shows that the fl = Jahn-Teller model is equivalent to an Ising model with 
interaction strength, 

J3,i = -9*Y, M *3,n^M l>n . (8) 
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The former argument is equivalent to the adiabatic elimination of bosons [5] . Adding 
the transverse field term, we find 

e- s H s e s = H s = ^ ]T (a+e' 2 £» ^^O^) + H.c.) , (9) 

3 

which shows that as long as g <C u) n , H s w H s . Previous works exploited this idea for 
the quantum simulation of Ising and Heisenberg models with trapped ions [I] . Here 
we are interested on the whole phase diagram of the model, where bosons suffer a 
displacement which cannot be neglected. 

We aim to give a mean-field description and to test its validity. Let us first quali- 
tatively discuss the problem. Mean-field theory may be well justified by two lines of 
reasoning that arise in two different fields: 

— Quantum magnetism.- The limits of large Q and large g should present a well- 
defined magnetic order in the x and z directions, respectively. Mean-field theory 
is a good approximation there. 

— Cooperativity effect.- In the limit of large t, we expect Q n to form a set of 
energies with a large energy separation Aui oc t. In the case that Aui 3> g, we may 
simply ignore high energy modes and keep n = only. We recover the celebrated 
Dicke model with infinite range interactions, for which mean-field theory is exact 
in the thermodynamic limit [3D]- However, as N — » oo we expect Au> — > so that 
this argument is only valid in a mesoscopic regime of spin-boson chains. 

The interplay between magnetic order and finite-size effects is a unique feature of the 
spin-boson model described by -ffjT- 



3 Mean-field theory 

Let us present our mean-field variational ansatz, which consists of a product state of 
spins in the x-z plane, and bosons displaced in the collective mode basis |21j . 

Wmy) =(g)|^-)®e-£"( a " a "- a " a ")|0) b , (10) 

3 

where \6j) — cos(0j/2)|O)j + sm(0j/2)|l)j. The mean-field energy is 
(^mfI^jtI^mf) = ^2 u n a* n a n -g^cos 9j (M* n a* n + M jin a n ) + — ^ sin 9j. (II) 

n j,n ' j 

We minimize the energy and get a set of coupled equations for the variational param- 
eters 9j and a n , 

(12) 

J U cos t = -q cot 0j, J ltj = 2 £ n 3? (f-M^M^) . 

Let us consider Periodic Boundary Conditions (PBC) to simplify the discussion, to- 
gether with the nearest- neighbours coupling ([3]). The problem is diagonalized by 
plane- wave modes such that 

uj n =wo + 2t(l-cos(2vrn/iV)), (13) 
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with wo the lowest mode energy; we have assumed constant local energies in @, 
ujj = ujq + 2t. The latter is a convenient parametrization so that the lowest collective 
energy mode is uiq, independent on the value of t. Collective mode amplitudes become 
simply 

M j<n = e~ iS Tr j /VN, (14) 

with n = 0, . . . , N — 1. It is worth noting that the eigenfunctions are normalized and 
that they follow the closure relation e z ~iv i ( m ^") = N5 m , n - 

By virtue of the PBC we can drop the index in 9j, 9i in the equations (|12[) . turning 
out that 



a n = — M*„ cos 9 — — cos 9— j= 



and that 



c 1 N ^ 

3 3 

^VN cos 9S n0 = ( " n *° = g °^ . (15) 



1 



-e N J — e 

N Q n 



E Ju = E ^ ( E M i^ M A = 2?* I e ^ 

l l \ n n / \n,l 

= 2 ^ (E i^^N8 n0 ) = ¥ := J- (16) 

Inserting (|16p in the mean-field Eqs. (|12[) we find the minimum at 

sin(0) = -n/(2J), \Q\/{2J) < 1, 

sin(0) = -1, |«|/(2J) > 1. (17) 

Thus, mean-field theory predicts a phase transition at a critical value |J? C | = 2 J . The 
behaviour of spins on both sides of the critical point can be defined by taking the two 
different asymptotic limits, 

cosff = 0-».|fi) = -^(|0)±|l)) ) 

(18) 

sin (9 = Q->\0) = {|0),±|1)}. 

The meaning of the two phases is obvious and corresponds to the discussion in the 
previous section. Note that the mean-field theory yields the same result that would 
be obtained by an effective Ising spin model. However, fluctuations around the mean- 
field solution are different in the original spin-boson lattice and thus, the validity of 
the mean-field description cannot be addressed with an effective Ising Hamiltonian. 




4 Gaussian fluctuations 

Quantum fluctuations are expected to diverge at the critical point. Near the latter our 
mean-field description is no longer valid. To quantify fluctuations we use a Gaussian 
approximation around the mean-field solution. As we are going to consider a particular 
case of the mean-field solution later on, which presents homogeneity in the spin states, 
we drop the index in 9j. 
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Let us define first fluctuations operators with respect to the bosonic degrees of 
freedom, Sa n — a n — a n . Spin fluctuations are more involved and we quantify them 
by using Holstein-Primakoff bosons [35]. We have to work in a spin rotated frame 
such that the mean-field state is the reference state for the Holstein-Primakoff trans- 
formation. To this aim we define the rotated operators a'j' y ' z , 



cos0ct| +sva6a], (19) 
cos 9a j — sin 9a j . 

The rotated Hamiltonian (from the original Hamiltonian ([2]) in the collective mode 
basis) becomes 




j n 

+ gJ2(cos6a] + sm9a?){M j<n a n + M* n al). (20) 

After the rotation, we can use the usual Holstein-Primakoff transformation [35] where 
the reference state is taken in the crj basis. The spin fluctuations are expressed by 

means of bosons bj , b\ , 



< aj = y/i - n sJ bj, (21) 
y t =2n Stj -l, 



where n s j — bjbj. It is straightforward to show that those operators indeed obey 
the commutation relations corresponding to the Pauli matrices. In the limit of small 
fluctuations we approximate 









(22) 


aj ~ bp 





The former approximations are only valid in the limit (bjbj) <C 1 which we shall verify 
later for self-consistency. Near the critical point, we expect that large spin fluctuations 
such that nonlinear terms in Eq. (|21j) become relevant. 

Substituting in the rotated Hamiltonian we arrive at the Gaussian fluctuations 
Hamiltonian, 

H G = Un&a) n 5a n + ^ Ajblbj 
n j 

+ Y J 9^0{MlJa\ l + M j<n 8a n ){b] + bj), (23) 
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where Aj := — fi sin (9 + 2g cos 9 J2 n (Mj. n ct n +c.c). Hg is diagonalized by means of 
a Bogoliubov transformation to spin-phonon fluctuation operators c m , 

H G = E ™cic m . (24) 

m=l,...,2N 

The new 2N bosons are related to the N boson and N spin fluctuation operators by 
a relation of the form 

5a n = ]T (W r ^ m c m + Vj%cl) , 

m=l,...,2JV 

b i= E (W^crn + V^cl), (25) 
m=l,...,2N 

where the matrices Wn%,, V„%, define a generalized Bogoliubov transfor- 

mation, which has to satisfy the condition that bosonic commutation relations are left 
invariant. We will find below a very convenient way to get an explicit result for that 
transformation. For the moment, we just show how to compute quantum fluctuations. 
Define first the vacuum of the eigenmodes of Ha, \Q) by the condition c m \ft) = 0. 
Then we define the variance per atom for a set of the original spin-phonon fluctuation 
modes and get the result 

F{Sa} = ^J2( n \ S ^ n \0) =£|V;W| 2 , 

n n,m 

F m = ^E<^W = El0 2 - ( 26 ) 

3 j,m 

In the next section we use these expressions to calculate quantum fluctuations in the 
case of translational invariant systems with PBC. 



5 Analysis of the problem under Periodic Boundary Conditions 

To impose PBC we substitute the eigenfunctions Mj jTl defined in Eq. (|14l) . We trans- 
form fluctuation operators to the plane- wave basis by the relations 5a j — ^ n Mj, n Sa n , 
bj = J2 n Mj.nbn, and use the relations (|I51 [TBI \T7} in the Hamiltonian (|23p . to get 

H G =^2uj n SalSa n + V g sin 6(5al n + 5a n )(bl + &_„) - V -^-xb\b n . (27) 
*■ — ' * — ' * — ' sin ft 

n n n 

Note that the closure relation in the plane- wave basis yields couplings between modes 
with opposite linear momentum, n and in terms of the form 8ci ri b— ni for example. 
To diagonalize He we work in the X, P representation. We define 



x ^ = vm( 6n -n + J Xbin = J^{bl n + K), 



(28) 



Note the relations Xl n — X^-n, Pj[ n = P^-n (/•* = a i b). Up to an additive constant 
corresponding to the vacuum energy, the problem can be stated as a sum of separable 
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Hamiltonians for each mode n, 



i i _ 2 i i / n \ 2 

Ha n — —Pan.Pa.—n + ~UJ „X n _ n X n _ — „ + — Ph.n.Ph. — n + 1 — I Xh.nX, 



- ffV^Tsm^j (X a ,„X b ,_„ + X b , n X a - n ) . (29) 
In matrix notation each of the Hamiltonians Hc, n can be written as 

H G , n = IJ2 Ki$X*nXv,- n + \ E P^nP^-n, (30) 



K (n) = 



2 



-1g\J fi(jJ n \ sin 0p 



where fj, and ^ run over labels a, &. 

Now we have to find a transformation defined by a matrix C/^™2 , X^ n = uj^}x^ n , 

x n,-n = J2 V ufyXv-n, P^n = J2 V Uf$Pv,n, P/i-n = T, u u t$Pv,-n that diagonal- 
izcs K {n \ and conserves the canonical commutation relations. The latter requirement 
is fulfilled automatically since the transformation is orthogonal due to the fact that 
j£{n) j g rea j anc j sy mm etric. 

The eigenvalues of K^ n ' can be split up into two branches, 



(^) 2 + ^±V 16 .9 2 ^l sin( 



| sine| 



2 N 2 



E±,n= ! 2 " ~- (31) 

The transformed eigenvectors X± >n and P± n are the columns of the matrix 

2 gv /ffi0„|sin8| E 2 _ ,„-(^|) 2 \ 

U [n) = I 2 °" 2 ; (32) 



there the normalization factor is 

= ^-,„ - (y^y) ) + 4 5 2 ^»l sin ^ ( 33 ) 

We write the new operators X± jn ,P± tn in second quantized form, 

X±,n = / I (4,-n + C ±,»)> 

V 2 / £^_ (34) 

1 /B±.» • 



From here we can easily get a transformation to relate <5a Si „, 5a„ to normal modes 
c± >n like in Eq. (|25|) and use those to get a final expression for the fluctuations defined 
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Fig. 1. Quantum Gaussian fluctuations for N = 20 ions. Black and blue lines correspond 
to phonons and spin-waves respectively. We have considered homogeneous chains and units 
such that Q = u>o = 1, t = 0.4. The critical point is at g c = 1/2. 



in Eq. @, 

F {Sa } 



1 1 



UJ,, 



n 

l l 



+ ,n 



P 2 n 2 , 

-,n \sin9\ 2 ) E- t 



1 



sin0| E-, n 
4 E+. n \sm9\ 



In the following we are going to analyze the results obtained by means of Eq. (|35|) . 
In Fig. [1] we plot the quantum fluctuations as a function of g and check that they 
diverge at the critical point. Although this divergence is a straightforward effect it 
will shed light on the scaling of fluctuations with N. We have chosen units such that 
the spin and boson fluctuations are of the same order. The mean-field critical point 
is defined by the condition fl = 2 J = Ag 2 /l)q. 

First we notice that E- ;n = if n = since E- Jn - — > c ■ n, with c a constant 
determining the speed of acoustic fluctuations at the critical point. Therefore from 
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(|35|) it can be shown that the the fluctuations are dominated by the term l/i?_ iTl . 
This fact, indeed, makes the fluctuations ill-defined for n = 0. On the other hand 
the contributions n ^ to 1/£L,„ in (|55|) are the source of a logarithmic divergence 
expected as N — > oo, which can be studied separately from the n = divergence, by 
looking at the well-defined modes with n ^ 0, 

1 A 1 1^1 I [ N dn f 1 dx , Ar 

sS^-ifS^»J~J?A ^*rC^r log,v - (36) 

We conclude that the amplitude of quantum fluctuations has to be understood by 
studying separately two different contributions: 

— Modes n — show the same behaviour as the one expected for the Dicke model: 
fluctuations diverge at the critical point, however, they are suppressed for g ^ 
g c . As N — > oo, the critical point becomes singular and any g ^ g c shows no 
fluctuations. 

— Modes « / contribute with a logarithmic divergence at the critical point due 
to the infrared divergence in the fluctuations. Out of the critical point, we expect 
an infrared cut-off in the fluctuation spectrum, and thus there is a well-defined 
N — > oo limit. 

Out of this discussion we get the following picture: Exactly at the critical point 
(g = g c ) fluctuations diverge and they are governed by the n = mode. Out of the 
critical point both n = and n^O modes contribute but the n = contribution will 
be suppressed in the thermodynamic N —¥ oo limit. 

Following those arguments, we calculate separately the total fluctuations (Fb, Fs a ), 
the contribution from n — 0, (-FL, Fsa a ), and the fluctuations from modes n ^ 0, 
(En/o\> J2n^o F$a n ) ■ In Fig. [2] (left), we plot the total fluctuations as function of 
g for different values of N. Note that, for concreteness, we consider the case Q = ojq, 
which leads to similar spin and boson quantum fluctuations. We find convergence in 
the large N limit. In Fig. [5] (right) we plot the n = contribution and check that 
it contributes to the divergence at the critical point, but it can be neglected in the 
noncritical region for large N, as expected from the single-mode Dicke model. In 
Fig. [3] we plot the n^O contribution and the total fluctuations as a function of the 
number of particles. We observe that fluctuations converge to the n^O contribution 
for large N, however, for the values used in this calculations there is an intermediate 
mesoscopic regime of N s» 10 particles, where fluctuations decrease with N. This 
effect is induced by the suppression of the n — contribution as N increases. Finally, 
Fig. @] shows the scaling of the n / contribution for different values of t. At the 
critical point (right), we check the logarithmic divergence estimated in Eq. (J3SJ) - Out 
of the critical phase, fluctuations converge to a steady value. A crucial observation is 
that for increasing t, the contribution from n^0 decreases, due to the energy cost 
of n ^ fluctuations. 



6 Conclusions and Outlook 

We have presented a mean-field theory for cooperative Jahn- Teller models that appear 
in a natural way in a variety of quantum optical setups. Our investigation has relied 
on the use of mean-field and spin- wave theory. The latter is limited in what concerns 
the description of the critical phase. However, we have shown that there is a regime 
of validity for mean-field theory that can be checked self-consistently by calculating 
the amplitude of quantum fluctuations. Our results show that the mean-field phase is 
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Fig. 2. Quantum Gaussian fluctuations for increasing number of ions N including all the 
Bogoliubov modes n (left) and considering just the n — mode (right). We choose values 
Qo = 1, fi = 1, t = 0.4. We depict only the spin-waves fluctuations ^ ? {5„}- 
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Fig. 3. Quantum Gaussian fluctuations for values Q = ujo — 1 and g = 0.6. We have chosen 
two values for the hopping: t = 0.4 (up) and t = 10 (down). 



determined by the interplay between the n = fluctuations typical from long-range 
Dicke models and the ra^O fluctuations that are well-defined in the thermodynamic 
limit. Our calculations show that quantum fluctuations decrease with N due to the 
suppression of the n = contribution and arrive at a steady value for large values of 
N. 

Our work is relevant to several experimental setups in quantum optics. In par- 
ticular, in the case of trapped ions, it was recently shown that a generalization of 
the cooperative Jahn- Teller model © could be implemented. This can be achieved 
either inducing interactions with lasers [1] or with magnetic field gradients [51] . Typ- 
ical energy scales are L)q, 0,g,t w 100 kHz, and ion chains of N = 2,... 50 ions. 
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Another promising system is circuit QED, where qubit-field couplings in the ultra- 
strong coupling correspond to values g rs u)q,Q in the GHz regime |7I31) . Recent 
theoretical proposal could allow one to induce sidebands in the qubit-field coupling 
to achieve a high degree of controllability of the parameters of the model [32]. In 
addition, our work could be complemented by the study of Jaynes-Cummings-likc 
qubit-field interactions, which also arise in a variety of quantum optical systems. 

An interesting issue that could be addressed by means of spin- wave theory is the 
effect of finite temperature. Experimental quantum optical systems like trapped ions 
and circuit QED systems are typically out-of-equilibrium systems. Ground states are 
created by adiabatic evolution, so that the correct physical description would corre- 
spond to an initial finite temperature state which evolves by adiabatically turning on 
some Hamiltonian parameter (for example, the spin-boson coupling, g). The theoret- 
ical description would thus rely on a non-equilibrium spin- wave theory to account for 
temporal evolution. Finite temperature effects would then lead to a finite number of 
gaussian excitations in the final state after the adiabatic evolution. 

We acknowledge QUITEMAD S2009-ESP-1594, FIS2009-10061, CAM-UCM/910758 
and RyC Contract Y200200074, and COST action "IOTA". 
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